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The inverse problem method is tested for a class of mean field statistical mechan- 



ics models representing a mixture of particles of different species. The robustness of 
the inversion is investigated for different values of the physical parameters, system 
sizes and independent samples. We show how to reconstruct the parameter values 
with a precision of a few percentages. 
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1 Introduction 

In recent years there has been an increasing interest in studying the inverse problem 
in statistical mechanics mostly due to the fact that the thermodynamic formalism on a 
macroscopic base has proved to be effective in a variety of scientific applications that span 
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from the investigation of real neural networks in biology [TTHT9] to behavioral ethology 
for flocks pp. In this paper we are interested in a particular class of models J3j[8] that have 
naturally emerged within the application of the statistical mechanics formalism to socio- 
economic sciences (see also [T(3ll2"U] and references therein). Their first and most elementary 
appearance can be traced back to the so called discrete choice theory proposed by Daniel 
Mc Fadden (15] after his celebrated success in predicting the number and distribution 
of customers of the Bay Area Rapid Transit before its construction. Discrete choice 
theory doesn't contain interaction between individuals and from the statistical mechanics 
point of view can be seen as a mixture of a finite number of discrete perfect gases; its 
inverse problem is mathematically elementary and its efficiency amounts to the proper 
identification of the different species of particles (see also [11]). The necessity to include 
the interaction among agents led W. Brock and S.N. Durlauf [2] to introduce, within the 
socio-economic context, the simplest interaction structure which is given by the mean field 
Curie- Weiss Hamiltonian model. In order to successfully generalize the discrete choice to 
the interacting case, it was defined in [3] a multi-species mean field model. In this paper we 
propose a robustness test of the inverse problem in the multi-species mean field case. We 
start from the knowledge of the exact solution of the model, both in the single populated 
system and in the bi-populated one, not only at the thermodynamic limit in analytic form 
but also at finite and increasing number of particles by accurate numerical approximations. 
This, together with the standard criterion of maximum likelihood, provides a relation 
between experimental and theoretical quantities and allows to tackle the computation of 
the free parameters of the model, namely interactions and magnetic fields, from observed 
data. 

We then generate the equilibrium configuration of the models, at different system sizes 
and for different values of the parameters i.e. interaction strength and magnetic field. By 
use of the inversion formulas we show how the reconstruction of the parameters is achieved 
and how his robustness depends on both system size and number of independent samples 
used. 

The paper is organized as follows. In section 2 we recall briefly the Curie- Weiss model 
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and we review how to solve the inverse problem in this single-population model. The 
generalization of such a model to systems composed of many interacting groups (the 
multi-species mean field model) and the solution of the corresponding inverse problem 
are presented in section 3. Section 4 presents and discusses a set of numerical tests 
for both the single and the bi-populated system for finite number of particles and finite 
number of samples. We first investigate how the average quantities, magnetization and 
susceptibility behave for increasing system sizes in the standard Curie Weiss model. We 
find, in particular, that while the magnetization is monotonically increasing, in agreement 
with the first Griffiths Kelly Sherman inequality [9j[10j[12] , the susceptibility has a mono- 
tonicity direction that changes with the values of the coupling constant with respect to 
its critical point. Both quantities reach their limiting value at the speed of the inverse 
volume. We then investigate how the experimental magnetization and susceptibility at 
fixed volume depend on the number of samples and stabilize when their number increases. 
The effectiveness of the inversion is tested for different values of the coupling constants 
and magnetic fields. The same procedure is applied to the bi-populated model and again 
the robustness of the inversion is tested for different values of the parameters. We find 
in all cases that the inverse method reconstructs, with a modest amount of samples, the 
values of the parameters with a precision of a few percentages. 

2 Inverse problem for the Curie- Weiss model 

Denoting with N the size of the population, the Curie- Weiss model is defined by the 
Hamiltonian: 



where <7j G {±1} is the spin of the particle i (individual), the parameter J > is the 
coupling constant and h is the value of the magnetic field. The joint probability of a 




(1) 
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configuration of spins cr = (<Ti, . . . , <7at) is given by the Boltzmann-Gibbs measure: 

exp(-H N (<r)) 



Pn,j,h{(t} 



exp(-H N (tr)) 



exp (^m 2 N (cr) + hm N (a] 

J2 exp ( N[±m 2 N {cr) + hm N (cr) 
ren N V v 



(2) 



where = { — 1, 1} and 



m N (o) = T7^^ 



i=l 

is the magnetization of the configuration cr. We point out that, in the inverse problem, 
the usual inverse temperature parameter is absorbed within the two free parameters J 
and h. 

We will denote by u(-) the expectation value with respect to Pn,j,}i- Heuristically, this 
distribution favors both the agreement of people's choices or opinions <jj, with an external 
influence h, and the agreement between individuals, being J positive (whereas for J < 
it would favor disagreement). 

The inverse problem amounts to compute the values of J and h starting from the 
knowledge of the magnetization average and variance. Of course, when dealing with real 
phenomenological data its solution is made in two steps. The first is the identification of 
the analytical inverse formula providing a possible explicit expression of the free param- 
eters (J and h) in terms of the mentioned macroscopic thermodynamic variables. The 
second is the evaluation with statistical methods of the macroscopic variables starting 
from real data. The problem is generically well posed because the unknown parameters 
(interaction and magnetic field) are as many as the measured phenomenological quantities 
(average magnetization and its fluctuation). 

Let start by observing that when h ^ and J > or h = and J < 1, the Curie- Weiss 
model satisfies the following property (see [1]) 

lim u(rriN(er)) = m(J, h) (3) 
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where m(J,h) is the stable solution of the model mean-field equation 



m(J, h) = tanh(Jm(J, h) + h). (4) 



By differentiating the identity ([3]) with respect to h we obtain: 

d 

i im ^T a; ( m ^( <7 )) = X 

w^oo ah 



where 



and 



dm(J,h) l — m 2 (J,h) 

X = dh = 1- J(l-m 2 (J,/i)) ^ 



d , , „ 9 /E ( r e n Jv m w( <T ) ex P(-^( cr )) 
— w(m^(cr)) = — 



<9/i E <Te n Jv ex P(-^( cr )) 

= jv(w(m^(<r))-w 2 (mtf(<r)))- (6) 

In particular, the right hand side of the last identity defines the finite size susceptibility 
Xn- By putting together (jSJ) and ([6]) we can compute the parameter J from the average 
value and the variance of the magnetization in the thermodynamic limit: 

j = 1 I (7 ) 

w 2 (mjv(r7)) ]hn N(u(m 2 N (tr)) - u*(m N (*))) 

N— >oo V / 

The external field h is obtained, in the large volume limit, by inverting the mean- field 
equation (BJ 

h = tanh -1 ( lim w(mjv(cr)) ) — J lim uim^icr)) (8) 

\N— >oo / N— >oo 

where J is given by (j7]). Formulas (jTJ) and (jSJ) solve the analytical inverse problem for the 
Curie- Weiss model as h ^ and J > or h = and J < 1. On the other hand, if h = 
and J > 1 the equation (jl]) has two different stable solutions ±m( J, 0) (see [4]). Thus, the 
identity fl3]) is not verified. In this case the inverse problem can be solved, by observing 
that there exists e > such that, whenever mjv(cr) G (±m(J, 0) — e, ±m(J, 0) + e), the 
following limit holds (see [5]) 

lim u(m N (cr)) = ±m(J, 0) (9) 

N— >oo 
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and then by applying to (Q the same procedure as shown above. The result is still given 
by formulas (J7J) and (jHJ) which conclude the analytical treatment of the inverse problem. 

For what it concerns the statistical part one has to provide an evaluation of the finite 
size average magnetization uj{m N (cr)) and susceptibility xn from the empirical data. 
We use the maximum likelihood estimation procedure. This method identifies the free 
parameters within a distribution by requiring that their value maximize the probability 
to obtain the given sample, under the condition that the sample is made of independent 
and identically distributed realizations of the random variables. 

Given a sample made of M independent spin configurations cr^\ . . . , cr( M ) all dis- 
tributed according to the Boltzmann-Gibbs measure (j2J), the maximum likelihood [71 [T5] 
function is defined by 

which, using the independence, can be rewritten as 

M 



i(^)=n%{ ffW } 



i=i 

M 

n 



To maximize the function L(J,h) we should compute the derivative of a product. Since 
a function and its logarithm reach the maximum in the same point, we consider the 
logarithm of the maximum likelihood function 

\nL(J,h) = Y / {-H N (tr®)-hi ^ exp(-H N (a-))j . 

The derivatives with respect to h and J of this function 



dkiL(J,h) 



dh 

i=i 



d\nL(J,h) _ N 
dJ ~ ~2 

i=i 



M, v 

i=i ^ ' 
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vanish for 



u(m N (cr)) 



u(m 2 N (cr)) 



1 M 



(10) 



M 



m 

i=i 



Therefore, the function L(J,h) reaches its maximum when the first and second momentum 
of the magnetization are calculated from the data according to ( jTOl) . The inverse problem 
is finally solved by the composition of ( TTOT) with (EJ) and flH]). In particular, denoting by 



and Xexp respectively the average magnetization and the susceptibility computed 



from the sample 



exp 



1 M 

i=l 



m N {cr 



Xexp 



m 



exp I i 



(11) 



i=l 



the estimators of the model's free parameters are 

1 1 



J, 



exp 



m 



exp 



Xexp 



h e xp — tanh (m exp ) 



J expfTlexp- 



(12) 
(13) 



3 Inverse problem for the multi-species model 

Since, within the applications we are interested in, we aim at generalizing the discrete 
choice theory [15] . we proceed toward the solution of the inverse problem for the multi- 
species mean field model. Formally such a model is an extension of the Curie- Weiss model 
to systems composed of many interacting groups. We consider a system of N particles 
that can be divided into k subsets Pi,...,Pk with Pi fl P s = 0, for I ^ s and sizes 
| Pi | = Ni, where ^4=1 Ni = N. Particles interact with each other and with an external 
field according to the mean field Hamiltonian: 

^ N N 

H n (<t) = -^Yl ■ i ':< C! ' a > ~ hi(Ti ■ ( 14 ) 

i,j=l i=l 
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The Oi G {±1} represents the spin of the particle %, while Jij is the parameter that tunes 
the mutual interaction between the particle i and the particle j and takes values according 
to the following symmetric matrix: 



No 



\ 





N 2 




Jll 


Jl2 




J lk 


Jl2 


J22 














J lk 


J 2k 




Jfcfc 



\ 



/ 



where each block J; s has constant elements Ji s . For I = s, 3u is a square matrix, whereas 
the matrix J^ s is rectangular. We assume Ju, J22, ■ ■ ■ , Jkk t° be positive, whereas Ji s with 
I 7^ s can be either positive or negative allowing for both ferromagnetic and antiferromag- 
netic interactions. The vector field takes also different values depending on the subset the 
particles belong to as specified by the following vector: 



iV 2 { 




where each h/ is a vector of constant elements hi. Indicating with mi(cr) the magnetization 
of the group Pi, and with «/ = Ni/N the relative size of the set Pi, we may easily express 
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the Hamiltonian (TEH) as: 



H n (ct) 



~ N {.2 aia s Ji s mi(cr)m s (<j) + ^aAmj(( 

l,s=l 1 = 1 



= -w(i(JD„m(o-),D a m(<r)) + (h,D„m(<r)) 

where m(cr) = (mi(<r),...,m fe (<r)), D a = diag{ai, . . . , a k }, h = (hi, 
the reduced interaction matrix 

' J\\ J\2 ■ ■ ■ Jlk 



(15) 



hk) and J is 



Ji2 J22 ■ ■ ■ J'. 



2k 



\Jlk J2k ■■• Jkk J 

The joint distribution of a spin configuration cr = (cr 1; . . . , un) is given by the Boltzmann- 
Gibbs measure Pn,3,\i related to the Hamiltonian f fl4|) . where again we consider the inverse 
temperature parameter (3 absorbed within the model parameters J and h. The well 
position of the model has been shown in [8]. In particular, in the thermodynamic limit 
the model is described by the following system of mean-field equations: 

k 



Jk 



m^J, h) = tanh ( ^ a iJu m i(3, h) + hi 
( k 

m 2 (J, h) = tanh I a iJil m i{3, h) + h 2 



(16) 



m k (J, h) = tanh f a iJik m(3, h) + h k \ . 

If the system (fl6|) admits a unique thermodynamically stable solution m(J,h) 
[mi(J, h), . . . , m fc (J, h)), the following identities hold (see [5]): 



lim uj(mi(a)) = m^(J, h) I = 1, . . . , k. 



(17) 



By differentiating the identities ffTTl) with respect to h s , s — 1, . . . , k, we obtain 



d 

J im ^7T a; ( m K cr )) = Xis l,s — 1, . . . ,k 
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where xis are the elements of the susceptibility matrix of the model. In particular, 

k 

dmi(3,h) d ( , /, >r-^ /T , N \\ 

X/s = — dh = dh V \ 1 + 2^ a P^P m p( J > h ) ) ) 

s p=i 

k 

= (1 - mf (J, h)) ^ s + at p Ji P x ps ) 

where <5; s denotes the delta of Dirac picked in I — s. Therefore, the susceptibility matrix 
X can be written as: 

X = P(I + JD Q x) (19) 

where P = diag{l — m1(J, h), . . . , 1 — ml(J, h)} and I is the identity matrix. Moreover, 
for each I, s = 1, . . . , k 

-u{mi(a-)) 



dh s y iy " dh s \ E^e-^W 

= A?" s ^u;(m;(cr)m s (cr)) - u{mi{a))u{m s {cr))\. (20) 



By computing the elements of x according to ( jl8l) and ( 1201) . by ([191 we get an expres- 
sion of the reduced interaction matrix J related to the average value and the correlations 
of the magnetizations in the thermodynamic limit: 

J = (P- 1 - X- 1 )^ 1 , (21) 

see [H]. Once the matrix J is determined, the elements of the vector h = (hi, . . . , hk) are 
obtained by inverting the mean field equations (JTBjl 

k 

hi = tanh -1 ( lim u(mi(cr)) ) — N ot s Jis lim u(m s (er)) I = 1, . . . , k. (22) 

s=l 

The previous formulas ( 12TT) and (1221) represent the analytical solution of the inverse 
problem when the system of mean-field equations (Tl6l) has a unique stable solution. When 
there are more stable solutions, identities (fTT|) have to be handled with care. Similar 
identities are, in fact, locally fulfilled around each solution and the solution of the inverse 
problem is still possible by applying to them the same procedure described above. The 
estimators are again given by (12T]) and ( 122]) . Before proceeding with the statistical part 
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of the inverse problem we notice that it is well posed since, in particular, the degrees of 
freedom of the problem are equal to k(k + 3)/2. 

Also in this case, starting from phenomenological data, we proceed with the help of 
the maximum likelihood principle. Consider a sample of M independent spin configu- 
rations cr^\ . . . , cr( M ) distributed according to the Boltzmann-Gibbs measure Pn,j,Ii] the 
maximum likelihood function related to the sample is 



L(J,h) = P N ,j, h {<T 



(1) „{M) 



M 

i=i 

M 



(0 



n 



exp(-H N ((T^)) 



i E (Te n JV ex p(-^( cr ))' 

Differentiating the logarithm of the likelihood function (|23|) 



lnL(J,h) = ^f -In £ exp(-^(<r)) 

1=1 ^ creUN 



with respect to ft./ and J; s , i, s = 1, . . . , k we obtain 



<91nL(J,h) A^a/Q; s 



i=l 



^ m((<r W )m> (i) ) - w(m i (cr)m s (cr)) 



These derivatives are equal to zero as the following equalities hold 
u(mi(cr)) 



1 M 

— E^(o-W) Z = 1, , As 



1 Af 

u{mi{a)m s {a)) = — J2 m ; (crW)m s (<xM) I, s = 1, 

M j = i 



fc. 



(23) 



(24) 



Therefore, the inverse problem for the multi-species model is solved by the composition 
of (|24|) with (l2Tj) and (|22|) . In particular, denoting by mi exp the average magnetization of 
each specie calculated from the data 

1 A/ 

mzezp = T7^2 m i(°' (t) ) ^ = 1, --- , A: 



i=l 
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and denned the matrices P exp = diag{l — ml exp , . . . , I— ml exp } and X exp , whose elements 
are 

( 1 M \ 

Xisexp = N s I — ^m/(crW)m s (cr w ) - m lexp m sexp 1 l,s = 1, . . . , k 

the model estimators are 

J exp — (P e xp — Xexp)^a (25) 

k 

ht exp = tanh _1 (m ; exp ) - ^ a sJu 

exp1Tl s eX p l = l,...,k. (26) 

8=1 

4 The inversion at finite volume and finite sample 
size 

When dealing with real data the elegant exactly solvable model has to be replaced 
by its finite size version. This is reflected both in the number of particles iV and in the 
number M of independent configurations in the sample, available from the statistical set. 
It is therefore important to see how the inversion formulas perform for different values 
of those quantities at assigned values of the parameters. The Curie Weiss model and its 
generalized multi-species version provide an ideal testing set not only because most of the 
applications concern mean field models but also because their finite size solution can still 
be handled thanks to the observation that the magnetization spectrum has a probability 
distribution that can be exactly computed. 

In this section we present a numerical test of our inversion procedure, both for the 
Curie- Weiss model and for its multi-species version. In both cases, for each choice of the 
size iV of system and of the free parameter values (J and h for the Curie- Weiss model, J 
and h for its multi-species generalization), the data that we are going to use are extracted 
from a virtually exact simulation of the equilibrium distribution. This is possible thanks 
to the mean-field nature of the models ([T]) and ffTBI) . which reduces the computation of 
corresponding equilibrium distribution to that of the weights of the 0(N) values of the 
magnetization. In this way, from Pn,j,h f° r the Curie- Weiss model and from Pjv,j,h for 
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its multi-species generalization, we can compute the finite size average magnetization and 
susceptibility and extract sequences of configurations. 

Although obvious, it is probably worth remarking that the parameter estimation 
method involves two approximations. The first one is in the inversion formulas (J7J),(JH]), 
( 12 ip and f )22|) . that require the infinite volume limit; the second one is the statistical 
error appearing in the evaluation of the averages and correlations through the maximum- 
likelihood estimators defined in ( fTOl) and ( 12^1) . In principle, the first approximation could 
reduce strongly the scope of the method to systems with very large number N of individu- 
als and it corresponds to estimating the finite size corrections of average magnetization and 
susceptibility. We don't go through this issue, rather we illustrate it with some numerical 
example to support the choice of the values of the parameters. Indeed, figure [TJ shows the 
finite size average magnetization u;(mjv(cr)) and susceptibility N{uj{m 2 N (cr))— u 2 {m^{(j))) 
for the Curie-Weiss model at different iV's for J = 0.6, h = 0.1 and for J = 1.2, h = 0.3, 
while the same quantities in the thermodynamic limit m(J,h) and x are represented by 
the horizontal lines (to ease the notation, in the figure and in the following we denote 
the finite size quantities respectively by m N and xn and we omit the dependences of the 
equilibrium magnetization). The figure highlights the monotonic behavior of and xn 
as function of N . In particular, is monotonic increasing for each value of the inter- 
acting parameter J, while xn is monotonic increasing as J < 1 and monotonic decreasing 
as J > 1. We point out that the different behavior of the finite size susceptibility is 
very useful dealing with empirical data because it tells us if the system is above or under 
the interacting parameter critical value before to apply the inversion procedure. Note 
that, for N > 5000 we have optimal approximations both for m and x i n the thermody- 
namic limit. The power-law fits in figure H] show evidence of the 0(./V _1 )-behavior of the 
finite size corrections both for magnetization and susceptibility, which entails the same 
0(iV _1 ) error in the estimation of J and h. Figure |3] gives evidence of the dependence 
of the estimators m exp and Xexp, given by (ITT)) , on the choice of the number M of the 
configurations of the sample that we use in the maximum-likelihood procedure. To asses 
the statistical dependence of m exp and Xexp on the sample {cr^\ . . . , cr( M )}, we computed 
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Figure 1: Finite size average magnetization tun (upper panels) and susceptibility xn 
(lower panels) as a function of N for the Curie- Weiss model for J = 0.6 and h = 0.1 (left 
panels) and for J = 1.2 and h = 0.3 (right panel). The blue continuous lines represent 
the magnetization m and the susceptibility \ in the thermodynamic limit. 

their values over a set of 20 independent instances of such samples. Thus from now on, 
and without ambiguity, we use the subscript exp to denote both estimators and their 
statistical mean over the 20 M-sample. We find numerical evidence that M > 10000 
stabilizes the estimations. In particular, the standard deviation of both m exp and Xexp 
as a function of M behaves as a power law M -0 5 (as J = 0.6 and h = 0.1, the fit of 
the standard deviation of m exp is aM~ a with a = 0.4933 ± 0.06, a = 0.013 ± 0.006, and 
goodness of fit R 2 = 0.9696, while those of the standard deviation of Xexp is bM'P with 
13 = 0.5175 ± 0.098, b = 3.269 G (0.7388, 5.8) and R 2 = 0.9343). 

In order to test numerically the inversion procedure for the Curie- Weiss model, we 
consider a sample of M = 20000 spin configurations {a^},i = 1, . . . , M, where tr® = 
(oi , ■ ■ ■ , crjy ) and iV = 10000. For a given couple of parameters (J, h), we extract the 
sample of M independent identically distributed spin configurations from the Boltzmann- 
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Figure 2: J = 1.2 and /i = 0.3. Upper panel: \m^—m\ as a function of iV together with the 
best fit aN b for the data in the right upper panel of figHJ We obtain a = 0.5047 ± 0.0037 
and b = —1.006 ± 0.002 with a goodness of fit R 2 = 1. Lower panel: \xn — x\ as a 
function of iV together with the best fit cN d for the data in the right lower panel of figiU 
We obtain c = 2.037 ± 0.019 and d = -1.006 ± 0.002 with a goodness of fit R 2 = 1. 

Gibbs probability distribution function Pn,j,h- Given (J,h), we consider 20 M-sample 
and we solve the maximum likelihood model for each M-sample independently; then we 
average the inferred values J exp and h exp of the model parameters, given by (1T21 and (TT3"|) . 
over the 20 M-samples. We consider J e [0.6,1.2] and h e [—0.3,0.3]. The obtained 
values for the case h = 0.1 and h = —0.1 are shown in fig JH and in figlH where J exp and 
h exp are plotted as functions of J. Note that the inferred values of the parameters are in 
optimal agreement with the exact values (continuous lines in fig3]and in fig ED- The J exp 
and h exp values generated for h = 0.3 are reported in figj6]to show that when the external 
field is enforced the prediction is quite good too: the points lie on the line of the exact 
agreement (red continuous lines in figED even if the error bars are bigger than in previous 
cases. 
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Figure 3: N = 10000, J = 0.6 and h = 0.1. Upper panel: Average magnetization 
m exp (blue crosses) as a function of M (number of the configurations in the sample) 
together with statistical error bars over 20 M-sample for the Curie- Weiss model. The 
blue continuous line represents the finite size magnetization for iV = 10000. Lower 
panel: Susceptibility Xex P (blue crosses) as a function of M (number of the configurations 
in the sample) together with statistical error bars over 20 M-sample. The blue continuous 
line represents the finite size susceptibility xn for N = 10000 
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h=0.1 




Figure 4: Left panel: J exp as a function of J for h — 0.1 (blue crosses). Error bars are 
standard deviations on 20 different M-samples of the same simulation (see text for the 
details of the simulation). The red continous line represents J exp = J. Right panel: The 
value of h exp (blue crosses) calculated from (JSJ) for the values of J exp in the left panel, as 
a function of J together with the statistical error over 20 M-samples. The horizontal line 
corresponds to the exact value of the magnetic field h — 0.1 
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As a test problem for the multi-species mean-field model we consider a system of 
N = 2000 particles divided into k — 2 equally populated subsets (Ni = N 2 = 1000) and a 
sample of M — 10000 independent spin configurations. Starting from 20 different couples 
of given values for the reduced interaction matrix 

Jn J12 



J \2 J22 

and for the external vector field 



Jn, J22 e [0.55, 1.2], J 12 G [-0.6, 1.1] (27) 



h u h 2 e [-0.3,0.3] (2* 



h 2/ 

we consider 20 M-samples for each couple (J, h) and we solve the maximum likelihood 
model for each one of them independently; then we average the inferred values J exp and 
h exp of the model parameters, given by ( |25i) and ( 1261) . over the 20 M-samples (as in the one 
population model). In figf7]the euclidean distances between 3 exp and the initial reduced 
interaction matrix J (blue stars) and between h exp and the initial external vector field h 
(red circles) are shown for each of the 20 choices of J and h (cases). We observe that, as 
in the one-population model (Curie- Weiss model), the inverse problem procedure to infer 
the experimental values for the coupling matrix and for the external field produces results 
in very good agreement with the initial values. In order to have a quantitative measure 
of the goodness of this procedure, we focus on the two most representative cases. 

If we consider the case 1, for which the distance between the initial data matrix J and 
the inferred matrix 3 exp gives the maximum value (see absolute errors in fig|7|), we have 

1.2 0.98\ / 1.173 ± 0.036 0.993 ± 0.028\ 

J exp = (29) 
^0.98 0.8 J yO.993 ± 0.028 0.794 ± 0.040^ 

and for the corresponding values of the external field 

/ 0.102 ±0.012\ 
K XP = \ • (30) 

\0.198± 0.011 J 

The errors on each value of the matrix J exp and the vector h exp are standard deviations 
across 20 different M-sample of the same (J, h)-simulation. 

18 




h=-0.1 




J J 

Figure 5: Left panel: J exp as a function of J for h = —0.1 (blue crosses). Error bars are 
standard deviations on 20 different M-samples of the same simulation (see text for the 
details of the simulation). The red continuous line represents J exp = J. Right panel: The 
value of h exp (blue crosses) calculated from (JS} for the values of J exp in the left panel, as 
a function of J together with the statistical error over 20 M-samples. The horizontal line 
corresponds to the exact value of the magnetic field h = —0.1 
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J J 



Figure 6: Left panel: J exp as a function of J for h = 0.3 (blue crosses). Error bars are 
standard deviations on 20 different M-samples of the same simulation (see text for the 
details of the simulation). The red continuous line represents J exp = J. Right panel: The 
value of h exp (blue crosses) calculated from (JSJ) for the values of J exp in the left panel, as 
a function of J together with the statistical error over 20 M-samples. The horizontal line 
corresponds to the exact value of the magnetic field h = 0.3 
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Figure 7: Upper panel: Absolute errors in reconstructing J and h. Distance between the 
reconstructed matrix 3 exp and the initial data matrix J (blue stars) and distance between 
h exp and the initial h (red circle) for 20 different choices of parameters J and h (cases). 
The values of 3 exp and h exp are averaged across 20 M-sample (see Section 4 in the text 
for the details of the simulations). Lower panel: Relative errors in reconstructing J and 
h. Maximum percentage errors for the reconstructed matrix J exp (blue stars) and vector 
h exp (red circle) for the same 20 multi-species cases considered in the upper panel. 
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The case 18 in figJTJ which gives the minimum value for the distance between J and 
3 exp corresponds to 

/ 0.601 ±0.022 -0.798 ±0.019 , 
J exp = (31) 
y-0.798 ± 0.019 0.901 ±0.020 

and to the external field 



, -0.201 ±0.005 , 
h exp = . (32) 

-0.300 ± 0.005 y 

The errors on each value of the matrix 3 exp and h exp are the standard deviations across 
20 different M-samples of the same (J, h)-simulation. 




5 Conclusions and perspectives 

In this paper we have tested the robustness of the inversion method in a class of 
statistical mechanics mean field models. The novelty of the results is both on the finite 
size behavior of the exact solutions and on the quality of the inversion for finite number 
of samples. Our findings show that with a modest investment on samples we are able 
to reconstruct the values of the parameters within a few percentages. The relevance of 
the problem comes from the necessity to have a fully tested method in the parameter 
evaluation from real data of socio-economic type, as started from the seminal work of 
Brock and Durlauf. The nature of the investigated model belongs to those without in- 
trinsic randomness but we plan to extend a similar analysis to those cases with random 
interactions like the Sherrington-Kirkpatrick model, and/or random network connections 
among agents, like in the diluted models. 
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